Imaging spin flows in semiconductors subject to electric, magnetic, and strain fields 
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Using scanning Kerr microscopy, we directly acquire two-dimensional images ol spin-polarized 
electrons flowing laterally in bulk epilayers of n:GaAs. Optical injection provides a local dc source 
of polarized electrons, whose subsequent drift and/or diffusion is controlled with electric, magnetic, 
and - in particular - strain flelds. Spin precession induced by controUed uniaxial stress along the 
(110) axes demonstrates the direct k- linear spin-orbit coupling of electron spin to the shear (off- 
diagonal) components of the strain tensor, txy. 

PACS numbers: 72.25.-b, 87.75.-d, 71.70.Ej, 72.25.Dc 



The abihty to control and measure electron spin de- 
grees of freedom in semiconductors has been proposed as 
the operating principle for a new generation of novel elec- 
trical devices with the potential to overcome the power 
consumption and speed Hmitations of conventional elec- 
tronic circuits Semiconductor devices utilizing 
electron spin generally require: i) transport of spin- 
polarized electrons from one location in the device to 
another, and ii) a means to manipulate the electron spin 
orientation, either directly with magnetic fields or indi- 
rectly with electric and/or strain fields that exploit the 
spin-orbit interaction by coupling to electron orbital mo- 
tion. Many spin-based semiconductor device proposals 
H, ^ H, 1^ are based on a field-effect transistor geom- 
etry in which electron transport occurs in essentially 2- 
dimensional structures. In order to design semiconductor 
structures whose function is based on electron spin it is 
necessary to understand the transport and flow of spin- 
polarized electrons, and how it is influenced by electric, 
magnetic and strain fields in these 2D structures. 

Using methods for scanning Kerr microscopy, we ac- 
quire 2D images of spin-polarized conduction electrons 
flowing laterally in bulk epilayers of n-type GaAs. The 
images directly reveal the spatial dependence of spin dif- 
fusion and spin drift in the presence of applied electric, 
magnetic, and - in particular - strain fields. Controlled 
uniaxial stress along the (110) axes induces spin preces- 
sion, revealing the direct (k-linear) spin-orbit coupling of 
electron spin to the off-diagonal components of the strain 
tensor e. The coupling may be characterized by an effec- 
tive strain- induced magnetic field B,;, which is shown to 
be orthogonal to the electron momentum k, and there- 
fore chiral for radially-diffusing spins. scales linearly 
with |k|, yielding a spatial precession of electron spins 
that is independent of electrical bias and is considerably 
more robust against the randomizing (ensemble dephas- 
ing) effects of spin diffusion as compared with precession 
induced by external magnetic fields. 

The samples are 1 /im thick, silicon-doped (n-type) 
GaAs epilayers grown by molecular beam epitaxy on 
[001]-oriented semi-insulating GaAs substrates. Doping 
densities are He— 1, 5, and lOxlO^Vcm^- 

rieces were 

cleaved along the [110] and [110] natural cleave direc- 
tions, and ohmic contacts allowed an in-plane, lateral 



electrical bias in the [110] direction. The samples were 
mounted in vacuum on the cold finger of an optical cryo- 
stat. To apply controlled uniaxial stress along the [110] 
or [110] axis at low temperature, the cold finger incorpo- 
rated a small cryogenic vise, whose lead screw was ad- 
justed via a retractable actuator. Since all the samples 
exhibited qualitatively similar strain- related effects, we 
show data from only the ne=10^^/cm'^ epilayer. 

A local, steady-state source of electrons, spin polarized 
along the sample normal ([001], or z), was provided by a 
circularly-polarized 1.58 eV laser focused to a 4 /xm spot 
on the epilayer. While this pump laser (typically 10-25 
/iW) injects 50% spin-polarized electrons and holes, the 
holes spin-relax and recombine rapidly, leaving behind 
a net spin polarization of the mobile conduction elec- 
trons 13] ■ These polarized electrons subsequently drift 
and/or diffuse laterally away from the point of genera- 
tion. 2D images of the resulting z-component of electron 
spin polarization were acquired by measuring the polar- 
ization (Kerr) rotation imparted on a linearly-polarized 
probe laser that was reflected from the epilayer surface 
and raster-scanned in the x-y epilayer plane. The probe 
beam (50-100 ^W) was derived from a narrowband tun- 
able cw Ti:sapphire laser, and also focused to a 4 /im 
spot. For lock-in measurement, the pump laser polariza- 
tion was modulated from left- to right-circular (injecting 
spins along ±z) at 51 kHz. 

The measured Kerr rotation is a strong function of 
probe laser energy near the GaAs bandgap (inset. Fig. 
la). Here, 30 /im separates the pump and probe spots, so 
that the signal arises solely from a non-zero spin polariza- 
tion of the electron Fermi sea. This energy-dependent re- 
sponse provides a relative (and in-situ) monitor of stress- 
induced bandedge shifts. When imaging, the probe en- 
ergy is tuned (as shown) to the broad maximum that 
exists below the 1.515 eV GaAs bandedge, to avoid per- 
turbation of the Fermi sea. 

Fig. la shows a 70x140 /im image of the steady-state 
electron spin polarization in the ne=10^^/cm'^ GaAs epi- 
layer due to spin diffusion alone. To emphasize smaller 
signals, the color scale has been adjusted so that white 
equals half the peak signal. The spatial extent of the 
measured spin polarization (~60 /zm edge to edge) is 
much larger than the focused pump laser (shown by the 
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FIG. 1: (a) 70 X 140 /im image of electron spin polarization in 
a 1 /im thick niGaAs epilayer (ne = 10^^cm~^) at 4K, acquired 
via Kerr-rotation (KR) microscopy. A circularly-polarized 
1.58 eV laser focused to a 4 /im spot provides a local, dc 
source of spin polarized electrons; this image indicates 2D 
spin diffusion. Inset: KR vs. probe energy near 1.515 eV 
GaAs bandedge. (b) With E=10 V/cm lateral electrical bias, 
showing spin diffusion and drift, (c) Cross-sections of spin 
flow vs. bias; dotted line shows 5.5 /im resolution. 



red spot), indicating diffusion of electron spins away from 
the point of generation. Images of radially-diffusing elec- 
trons provide information on spin flows along all k di- 
rections in the x-y sample plane, which will prove useful 
later in confirming the direction of B^. The characteris- 
tic diffusion length is given by the electron spin lifetime 
Tg (measured independently via ultrafast techniques) and 
spin diffusion constant D^. Fits to a 2D drift-diffusion 
model (described later) indicate Dg ~3 and 15 cm^/s 
for the ne=l and 5xl0^^/cm'^ samples at 4K respec- 
tively, values in accord with the charge diffusion con- 
stants. Spin-polarized electrons can also be induced to 
drift laterally along an applied electric field E (Fig. lb). 
This "spin drag" , as revealed by Kikkawa and Awschalom 
0, occurs over length scales >100 /im. 2D images of spin 
drift and diffusion (Fig. lb, where E=10 V/cm) show 
roughly elliptical contours of constant spin polarization, 
with major axis determined by electrical bias and lat- 
eral spin flow extending beyond 150 /xm. Fig. Ic shows 
normalized line cuts through a series of images. 

Spin-orbit coupling in GaAs permits coupling to elec- 
tron spin degrees of freedom through the spatial part of 
the electron wavefunction, thus allowing induced preces- 
sion of electron spins without external magnetic fields 
d, 0, 01 . Spin-orbit effects lead to spin splittings of 
the conduction band along particular crystal momenta 
k, and can be characterized by effective magnetic fields. 
Bulk inversion asymmetry (BIA) arises from the lack of 
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FIG. 2: (a-c) 80x80 fim images of 2D spin flow (E=10 V/cm) 
at 4K, showing induced spin precession with increasing [110] 
uniaxial stress, (d) KR vs. probe energy for the images, 
showing blueshift of GaAs bandedge with stress, (e) Line 
cuts through the images. Inset: Spatial frequency of spin 
precession vs. bandedge shift. 



inversion symmetry in GaAs, leading to a spin-splitting 
for electrons with k along the (110) axes (but no split- 
ting along (111) or (100) axes). This ubiquitous coupling, 
cubic in |k|, is the origin of the D'yakonov-Perel' mecha- 
nism of electron spin relaxation in bulk GaAs. A second 
spin-orbit effect arises from structural inversion asymme- 
try (the SIA or "Rashba" term |0), as typically found 
in 2D heterostructures. Inversion asymmetry of the con- 
fining potential along the growth direction z (typically 
[001]) can often be characterized by an electric field E^, 
giving a Rashba Hamiltonian Hu cx cr -(k x E^). For 
typical 2D heterostructures with k in the x-y plane, the 
effective magnetic field k x E^ is therefore in-plane and 
orthogonal to k, with magnitude linear in |k|. Control of 
the Rashba spin-orbit term, through a gate-tunable E,, 
is the basis of the original Datta-Das spin transistor ^ . 

Figure 2 demonstrates an additional spin-orbit effect; 
namely the coupling of the electron spin to the strain 
tensor e. As described previously 0,0^3)^3) stress 
along the (110) axes of GaAs induces k-linear spin split- 
tings in the conduction band through the off-diagonal 
(shear) elements of e. The strain Hamiltonian is Hs = 
C3(T • tp, where {ip^,ipy,(p^) = {txyky - exzkz.f-yzkz - 
^yxkx, ^zxkx-ezyky), {x, y, z) are the principle (100) crys- 
tal axes, and the constant C3 depends on the interband 
deformation potentials. Stress applied along the [110] or 
[110] axis of GaAs gives in-plane shear Cxy = Cyx 7^ . 
Thus for electrons moving in the x-y plane, Hs has simi- 
lar symmetry to as it describes an in-plane effective 
magnetic field, orthogonal to k, with magnitude linear in 
|k|. 

Figs. 2a-c show 80x80 /im images of steady-state elec- 
tron spin flow (k II [110]) in the presence of increasing 
[110] uniaxial stress. Spin precession is observed, with 





FIG. 3: (a-c) 80x80 fim images of 2D spin flow (E=10 V/cm) 
at 4K, with increasing applied magnetic field Bapp=3.5, 9, and 
32 G along [110]. (d) Line-cuts through the images. Inset: 
Comparing line-cuts through Figs. 2c (black) and 3c (red). 



increasing spatial frequency, indicating a strain-induced 
effective magnetic field oc ea;j,|k|, oriented along [110]. 
e-xy is inferred from the measured blueshifts (see Fig. 2d) 
of '^1 and 2 meV in Figs. 2b and 2c respectively, indi- 
cating strain \exy \ ~1.5 and 3.0x10"^ (and applied stress 
~3.6 and 7.2x10* dynes/cm~^ These strains are 

small compared to typical ^1% strains due to lattice mis- 
matched growth; in fact, considerable care was required 
in sample mounting to avoid spurious and inhomogeneous 
strains during cool-down. With the cryogenic vise, the 
stress-induced precession of electron spins is controllable, 
reversible, and uniform over the sample. Line cuts along 
[110] (Fig. 2e) show many precession cycles (> Stt ro- 
tation). The inset of Fig. 2e confirms that the spatial 
frequency of the induced precession (oc Bj) scales linearly 
with the observed bandshift (oc ^xy)- 

Only shear (off-diagonal) strain gives k-linear spin- 
orbit coupling to electron spins. Strain along the princi- 
ple axes, either applied or arising from, e.g., lattice mis- 
matched growth along [001], should not influence elec- 
tron spins to leading order. However, shear strain and k- 
lincar coupling should exist in lattice-mismatched struc- 
tures grown along [110] or [111]. Because shear strain re- 
duces the symmetry of the zincblende GaAs crystal, these 
strain effects can be thought of as arising from BIA in a 
crystal of lower symmetry. In lower symmetry wurtzite 
crystals such as CdSe or ZnO, BIA alone gives k-linear 
spin splitting of electrons One should not think of 
Be as arising from an electric field {e.g., a stress-induced 
piezoelectric field ,17]), since any such field is screened in 
bulk metallic samples. 

The line-cuts in Fig. 2e show that dc spin flows, pre- 
cessing due to strain, are in phase at large distances 
from the point of generation. This robust behavior is 
in marked contrast with the rapid spatial dephasing that 
occurs when external magnetic fields are used to induce 
precession in dc spin fiows, as shown in Fig. 3 Figs. 3 
a-c show spin flows in the presence of an increasing ap- 




FIG. 4: 50x50 fim images of 2D spin diffusion (E=0) at 4K. 
(a) Stress=0, Bapp=0. (b) Stress=0, Bapp=16 G oriented 
along [110] as shown, (c) Bapp = 16 G, with [110] uniaxial 
stress. Spins diffusing to the right precess, while those diffus- 
ing to the left do not. Thus Be is chiral for radially-diffusing 
electrons (see diagram), (d-g) Maintaining [110] stress, B^pp 
is rotated by 180 degrees in-plane, (h) Stress is switched to 
[110], reversing chirality of Be. 



plied magnetic field (Bapp \\ [110]), with line-cuts in Fig. 
3d. As Bapp increases, the ensemble spin polarization be- 
comes dephased at distances beyond one precession pe- 
riod, particularly when the precession period falls below 
the spin diffusion length. This pronounced spatial de- 
phasing of dc spin flows is due to the randomizing nature 
of diffusion. The net spin at a remote location is the 
combined sum of many random walks. Each path takes 
a different amount of time, giving a different degree of 
spin precession, leading to dephasing. Future spintronic 
devices based on magnetic fleld manipulation of diffusive 
spin flow may thus be practically limited to a regime re- 
quiring TT rotation or less. 

A comparison of linecuts through Figs. 2c and 3c are 
shown in the inset of Fig. 3d. Clearly, the spatial coher- 
ence of dc spin flows is preserved over more precession 
cycles (and greater distance) when the spins are manip- 
ulated with strain instead of magnetic fleld. This is a di- 
rect consequence of the ]k]-linear nature of B^, which cor- 
relates precession frequency with electron velocity (and 
therefore position). Indeed, if electrons moved along only 
one dimension, the spin flow would not dephase at all (it 
would still, of course, decohere). Scattering from -|-k to 
-k would simply reverse the direction of precession, lead- 
ing to an exact correspondence between spatial position 
and spin orientation. 

The images in Fig. 4 conflrm that B^ is orthogonal 
to k, such that B^ circulates around the point of in- 
jection for radially diffusing electrons. Image (a) shows 
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FIG. 5: 50x80 ^im images of 2D spin fiow at 4K with in- 
creasing |k| (current), (a-c) With [110] stress and B^pp^O. 
Spatial period of precession is fixed, (d-f) Stress=0, Bapp=6 
G; spatial period varies. Images c') and f) are simulations of 
c) and f), using the model described in text. 

unperturbed spin diffusion (no stress and 'Bapp—0). In 
(b), Bapp=16 G along [110] as indicated. Spins precess 
uniformly, regardless of k, giving a faint annulus of op- 
positely oriented spins (negative signal) surrounding the 
injection point. In (c), stress is applied along [110], and 
the image becomes asymmetric. Electrons diffusing to 
the right along k || [110] undergo precession, while those 
diffusing to the left along -k do not. That is, the to- 
tal field (BT=Be-|-Bapp) is finite for spins diffusing to 
the right, but Bt is effectively zero for spins diffusing to 
the left. This image is consistent with a uniform B^pp 
added to a circulating B^, as shown. Maintaining [110] 
stress, the image asymmetry rotates and ultimately re- 
verses as Bapp is rotated in the x-y plane through 180 
degrees (c-g). Finally (h), when stress is applied along 
the [110] axis, the asymmetry again reverses, indicating 
the opposite chirality of B^ . 



For comparison with data, we derive and solve the spin 
drift-diffusion equations that describe dc spin flows in the 
presence of electric, magnetic and strain fields. For sim- 
plicity, the [110] strain axis is taken here to be the x-axis, 
and the electric and magnetic fields are in the x-y sample 
plane. The spin polarization is described by the (px, Py, 
Pz) components of a 2x2 density matrix, where pz gives 
the ensemble spin density. The equations are Oipx = 
-O2PZ, Oipy ^ -Ozpz, and O^pz-Oipx-OzPy = -Gz, 
with operators d = D\/^+p'E-V-(Cse)^D-l/T2, O2 = 
-CBBy + Cse{2DVx+f^Ex), O3 = CBBx + Cse{2DVy + 
liEy), and 04 = + pE-V- 2{CsefD - 1/Ti. D is 
the spin diffusion constant, p, is the electron mobility, T2 
and Ti are transverse and longitudinal spin lifetimes, the 
magnetic coupling constant is Cb = the strain cou- 
pling constant Cs is given in Ref. (15i] , Gz is the spin gen- 
eration rate, and e is the off-diagonal strain element t^y 
Drift and diffusion appear in the strain coupling terms, 
due to its dependence on k. The 1 pva thick n:GaAs layer 
is small compared with a lateral spin diffusion length, 
thus, the problem is 2-dimensional. The equations are 
solved with numerical Fourier transform methods. 

In contrast with the case of an applied magnetic field, 
both the model and the data reveal a consequence of k- 
linear B^: the spatial period of precession is independent 
of applied electrical bias. Fig. 5 confirms this, show- 
ing spin flow along [110]. In the presence of strain (Fig. 
5a-c) , the spatial precession period is independent of elec- 
trical bias, whereas for the case of applied magnetic field 
(d-f), the spatial precession period clearly increases with 
increasing electrical bias. Functional devices based on 
rotation of spin from one point in space to another {e.g.^ 
from source to drain contacts in spin transistor designs 
0,0,13, 6]) may well benefit from the freedom to operate 
at variable electrical bias. Along with spin manipulation 
via Rashba coupling (also linear in |k|), the spin-orbit 
coupling of spin flows to shear strain, as detailed in this 
work, also affords this flexibility. 
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